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1. BACKGROUND 


1.1 Problem Description 


Energy costs have risen dramatically in the past several years. 
This has impacted all segments of aviation. Operations that 
have historically reflected time costs, must now be reoriented 
toward minimization of a cost function heavily dominated by fuel 
costs. Many accepted flight planning practices are being 
reevaluated and redesigned. 

One useful tool in this effort to reflect the realities of fuel 
costs is a model to estimate fuel burn. 

1.2 Fuel Model History 


MITRE has developed for the Federal Aviation Administration 
(FAA) an approach to fuel burn modeling (References 1, 2). This 
model estimates fuel consumption as a function of path profile 
(geometric) data. The Internal aircraft parameters such as 
thrust are not required for the model to estimate fuel 
consumption. 

The resultant FAA Fuel Burn Model is a flexible, efficient tool 
for many applications. It can be used to evaluate proposed 
alternative routing structures, and for many other tasks 
(References 2, 3). It also has the capability to be adapted to 
a flight planning role. 

1.3 Project Concept 


MITRE was asked by the FAA and NASA’s Langley Research Center to 
develop a concept for a General Aviation flight planning aid. 
This aid was to help provide the general aviation (or commuter) 
pilot with a first cut, fuel efficient, flight planning aid. 

The aid had to be relatively Inexpensive in concept, and easy to 
use, to fit the general aviation budget and operating 
environment . 

The basic concept revolved around the use of the fuel burn model 
on a hand held programmable calculator. This would permit the 
pilot to use it in an airborne update environment, as well as in 
a ground planning environment. In the ground environment, the 
calculator might be linked to other computational devices, to 
identify weather data, or for initialization of the data base 
for in-flight use. 


1«4 Solution Concept 


The solution approach was to take the existing fuel burn model 
as a basis for calculation of fuel burns over specified 
segments. The optimal path would be established by a dynamic 
programming technique (see for example Reference 4). This 
technique iteratively, but systematically, constructs the set of 
feasible flight paths and iterates to the single optimal. 

1.5 Report Organization 


Section 2 will review the salient facts on the basic fuel burn 
equation, as applied to general aviation aircraft. Detailed 
derivations are contained in Appendix A and B. Section 3 
describes generlcally the input data necessary to the flight 
plan optimization. A more detailed description of a possible 
operational scenario is found in Appendix C. Section 4 presents 
the basic dynamic programming approach. Section 5 then reviews 
the interaction with the pilot, in terms of implementation of 
the optimal plan, and in terms of the process of in-flight 
update. Finally, the specialization of the dynamic programming 
concept for a hand held programmable calculator is addressed in 
Section 6. 
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2. DEFINITIONS AND ASSUMPTIONS 


2.1 Assumptions 


Several assumptions inherent in the development of the fuel burn 
model are continued here in the application to fuel efficient 
flight paths. The primary of these are: 

o Weight changes over Increments will not materially 
affect fuel burn calculations. 

o Flight path angles are small: cos (Flight path angle) 
approximates 1. 

o Velocity and altitude changes within an increment are 
linear. 

2.2 Fuel Burn Equation 


The fuel burn equation has been developed to measure fuel 
consumption as a function of path profile characteristics. The 
equation accepts as Inputs (see Table 2-1): 

o Stored constants and functions 

- general non-specific functions (density, etc.) 

- specific constants related to aircraft and engine 
type 

o Path/profile increment data 

- altitude 

- velocity 

- time or distance 

- weight 

- configuration (flaps, gear) 
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TABLE 2-1 

FUEL CONSUMPTION ESTIMATION ALGORITHM 
INPUTS/OUTPUTS 


Inputs 

• Stored Constants 

— General Nonspecific Constants (Density, Etc.) 

— Specific Constants Related To Aircraft and Engine Type 

• Path/Profile Increment Data 

— Altitude 
— Velocity 
— Time 
— Weight 

— Configuration (Flaps, Etc.) 

Outputs 

• Calculated Fuel Burn Over The Increment 

• Performance Limit Checks 



The fuel burn equation has been derived from basic principles of 
conservation of energy. Energy balance requires: 

(Energy In) - (Energy Loss) 

= (Energy Change) 

Specializing this to aircraft flight, we have: 
f (Thrust) “ f(Drag) 

- f (Kinetic Energy Change) 

+ f (Potential Energy Change). 

The details of the derivation are presented in Appendix A. 
Additional material on jet aircraft derivations is found in 
Reference 1. The initial derivation is given in Figure 2-1 for 
ISA conditions. Figure 2-2 gives the construction of the 
R 2 multipliers reflecting flaps and gear in low altitude 
flight. Figure 2-3 shows the modified equation necessary to 
reflect actual wind data. It requires mainly a correct 
accounting of the two factors: 

o Aircraft flight through the airmass (leading to 
instantaneous fuel flow rate) 

o Aircraft progress along the ground track (leading to 
real forward progress along trip itinerary) 

For simplification of the following discussions, we will simply 
talk of fuel burn (B) as calculable as a function of W, h, V, 
and T with aircraft constants being known and included. The 
effects of actual winds, flaps, and gear will be noted only 
when they have a significant impact on the fuel burn. 
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FIGURE 2-1 

FUEL BURN EQUATION FOR TURBOPROP (ISA CONDITIONS) 
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FIGURE 2-3 

FUEL BURN EQUATION FOR TURBOPROP 
(ACTUAL WINDS AND TEMPERATURE) 



3. INPUT DATA 


The general solution approach will be developed In Section 4. 
This section first defines the content and source of the data 
that is needed as input to the dynamic program. 

3.1 Stored Data 


Stored data will be specific to the individual aircraft type. 
Data will Include: 

o Aire raft /engine type - (to identify data) 
o Aircraft specific constants for the fuel burn equation 

o Maximum service ceiling (feet) 

o Maximum velocity (knots) 

o Maximum takeoff weight (pounds) 

o Operating empty weight (pounds) 

In addition, stored data for the Individual aircraft or operator 
could include : 

o Effective service ceiling 

o Adaptations to the aircraft type generic constants 
3.2 Pilot Input Data 


The flight specific data is supplied by the pilot. This would 
normally include the following: 

o Aire raft /Engine ID (to verify stored data Is 

applicable) 

o Departure and Arrival Velocities (knots) 

o Departure and Arrival Altitudes (feet of density 

altitude) 

o Effective Service Ceiling (feet) (for this trip, 

reflecting availability of equipment, known operator, 
or ATC imposed constraints) 

o Landing Gross Weight (pounds) (defined as empty weight 
plus payload and reserve fuel and reflecting other 
operator imposed constraints) 
o Number of Waypoints (generally keyed to ATC defined 
navigational aids, fixes or airports) 

For each of the trip waypoints, the pilot (or an automated link 
to a weather data base) must provide: 
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o Course Heading (degrees) 
o Magnetic Variation (degrees) 

o Distance to Waypoint (or X, Y Ground Coordinates) (nmi) 
o Number of Altitudes (for which weather data is to be 
input) 

For each waypoint, and for the selected number of altitudes, the 
following are input: 

o Density Altitude (feet) 
o True Wind Direction (degrees) 
o Wind Velocity (knots) 
o Temperature (^F) 

3«3 Processing of Weather Data 


The weather data as input must be processed to a form more 
suitable to the dynamic program. The following conversions are 
made : 


o Wind Direction from True to Magnetic: 

Wind (Magnetic) 

= Wind (True) + Magnetic Variation 
- Wind into North-East Components : 

WN = Wind Velocity [Cos(Wind Magnetic)] 
WE = Wind Velocity [Sin(Wind Magnetic)] 

where 


WN = North Wind Component 
WE = East Wind Component 

Several interpolations are determined. Winds and temperatures 
at a waypoint are fit over the several altitudes by a least 
squares regression: 
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WN^ = CjL + C2h 
- C3 + c^h 
Temp^ = C 3 + C^h 



where 


WN|^ - North wind component at altitude h 
WEj^ = East wind component at altitude h 
h = Density altitude 

Temph = Temperature at altitude h (°F) 

= Constants from regression 

These values then may be linearly interpolated as necessary 
between the pilot input defined waypoints. Once the necessary 
interpolations are completed, the winds at each node point are 
converted to head wind components as follows: 

WCh = C7 + Cgh 

where 

C7 = C-|^ cos (course) + C3 sin (course) 

Cs = C2 cos (course) + C4 sin (course) 

WC^ = headwind at altitude h 
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4. MATHEMATICS OF DYNAMIC PROGBIAMMING MATRIX 


This section presents the basic dynamic programming approach to 
the derivation of a fuel efficient profile. This approach uses 
the ground track and wind /temperature data described in Section 3. 
It presumes the filed ground course is followed. The resultant 
is a fuel efficient altitude and velocity at each user defined 
waypoint . 

4.1 Overview 

For the purpose of illustrating the approach, we will assume for 
the moment that velocity is known at each waypoint. The basic 
dynamic program approach is illustrated for a simple four-step 
problem in Figures 4-1 through 4-6. 

We initially define a set of potential density altitudes for 
each of the defined distance waypoints (Figure 4-1). The 
altitudes for the departure and arrival points are known, and 
thus some immediate simplification is possible (Figure 4-2). We 
start from the arrival point and iteratively, but selectively, 
calculate fuel burns. First we calculate the fuel burn from 
each candidate altitude at distance d 3 to the arrival point 
^4 (Figure 4-3). We then consider distance d 2 » For each 
candidate altitude at d 2 > we calculate the fuel burn to each 
altitude at d 3 . Since we already know the fuel burn from d 3 
to d^, we can calculate the most fuel efficient route from 
each altitude point at d 2 to d 4 . The solution for one d 2 
altitude is shown in Figure 4-4. Iterating this for each 
altitude at d 2 we have the most fuel efficient path for each, 
as shown in Figure 4-5. Note that one of the altitudes at d 3 
does not lie on any fuel efficient (partial) path, and thus is 
discarded from further consideration. We finally calculate the 
fuel burn from the single altitude at the depart point d^ to 
each altitude at d 2 * We already know each (partial) fuel 
efficient path from altitudes at d 2 to d^, so we now can 
calculate the optimal from d;^ ^4» shown by the heavy line 

in Figure 4-6. 

4.2 Definitions 


In f.act, we need to determine both altitude and velocity at each 
distance point. The same dynamic program approach holds, as 
shown in Figure 4-7. The primary change is in the number of 
fuel burn calculations that are required at each step. Thus for 
the example shown, at d 3 , we have 4 altitudes x 5 velocities = 

20 values (versus just 4 altitudes). 
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• Define Distance and Aititude Points 

— Limited in Number by Computational Capacity (4 Is Example) 
— May Correspond (in Part) to User Defined Waypoints 
— Divisions Need Not Be Linearly Chosen 

• (Velocity is Treated Analogously) 


t • • • • 

Altitude 

I • • • • 


Depart 0 



Distance 




0 Arrive 


FIGURE 4-1 
MATRIX DEFINITION 



• Altitude at Depart/Arrive Points Is Known 


X 

• 

• 

X 

X 

• 

• 

X 

X 

• 

• 

X 

Depart % 

# 

• 

Arrive 


FIGURE 4-2 

MATRIX SIMPLIFICATION 



• Calculate Fuel Burn 

— From Each of Candidate Altitudes at Distance Point d3 
— To Arrive Point 

(Assuming for Moment that Velocity is Known) 


Depart % 
di 




d2 



Ln 


FIGURE 4-3 

DYNAMIC PROGRAM-ds STEP 



Calculate Fuel Burn 

— From Each Candidate Altitude at Distance Point d 

— To Each Candidate Altitude at Distance Point d3 

Calculate 

— For Each Altitude at 62 

— Smallest Sum of Fuel Burn d2 to d3 to d4 

Save Only Minimum Fuel Burn Path 

— For Each Altitude at d2 



FIGURE 4-4 

DYNAMIC PROGRAM-dj STEP 


For Each Altitude at 62 , Have Fuel Optimal Path to d4 

Note Some d3 to d4 Path Segments Deleted From Further 
Consideration 



FIGURE 4-5 

DYNAMIC PROGRAM-dj STEP RESULT 



Repeat Previous Step for d-j 

Note Only Single Candidate Altitude at d 

Complete Optimal Path Is Defined 



FIGURE 4-6 

DYNAMIC PROGRAM-di STEP 





We now proceed to introduce the basic notation we will need. We 
define a three dimensional matrix in 

o Density altitude (index i) 

o Velocity (index j) 

o Distance (index k) 

We define the limits of (i, j,k) as M, N, P respectively, not 
necessarily equal. 

We will make the definition of altitude and velocity with ranges 
independent of distance. (This assumption can easily be 
changed, as noted later.) 

Altitude elements are defined by the user or as 

“i = Hmin + SM)<Hmax - 

where 


^MIN " minimum altitude 

(departure or arrival) 

^MAX ~ effective service ceiling 

(if less than maximum service celling; otherwise 
maximum service celling.) 

Similarly, velocity elements are defined by the user or as 
“ %IN 

where 


^MIN ~ minimum velocity 

(departure or arrival) 

VjjE ““ velocity not to be exceeded 


We will define distance elements by the user waypoints. If the 
number of waypoints is less than P (limit of distance index) 
then additional distance elements will be added. This will be 
accomplished by adding a distance element D^ such that 


D 


k 



+ D 

q 


where D is chosen 

q 
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such that Dq+i ~ t)q = Max„ (Dm+1 " 

At this new distance element the weather data is linearly 
Interpolated between distance elements and and the 

course heading is set equal to the course heading at Dg. This 
is repeated until the number of distance elements equals P- 

Finally, the arrival weight is defined as the weight at point 
(1,1, P) in the matrix. This weight must be greater than the 
specified aircraft operating weight empty. 

4.3 Matrix Steps 


The number of calculations implied by the matrices defined as in 
Figure 4-7 will expand rapidly for relatively small numbers of 
distance, altitude, and velocity points. Thus it becomes 
desirable to define the points carefully so they have relevance 
to the final optimum flight path. It is also important to limit 
calculations when they are clearly not relevant to the final 
path. Table 4-1 lists some of these considerations. The limits 
of aircraft performance can be used to reject some paths. 
Similarly smoothness rules can be applied. We know, for 
Instance, that the optimal path will not have a sequence of 
significant climb and descent combinations. These 
considerations are used in the logical definition of matrix 
steps below. 

1. Initial Segment 


a. Perform calculations between distance points Dp to 
Dp_i. Calculate fuel flow rate from each 
(i,j,P-l) to (1,1,P). Insure this is no less than 
idle fuel flow. 

b. Set path flags to indicate path to the arrival point: 

OH(i,j,P-l) = 1 
OV(i.j.P-l) = 1 

c- Set flags to indicate aircraft is descending towards 
arrival point and decelerating: 

ALT(i,j,P-l) = 1 

VEL(i,J,P-l) = 1 
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TABLE 4-1 

SIMPLIFICATION OF CALCULATIONS 
(RESTRICTION OF SEARCH) 


• Performance Constraints 

— Maximum Altitude 
— Maximum Velocity 
— Maximum Rate of Climb 
— Maximum Rate of Descent 

• Smoothness Rules 

— Climb After Descent 
— Velocity Increase After Decrease 

(Limits Imposed by Wind and Temperature) 
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These flags will be set to 0 when aircraft is 
ascending from departure point and accelerating 
respectively. 


d. For the next step, set 


\ " °p-i* 
2. Generic Segment 


Perform calculations between two distance points 
®k-l to Dfe. 

The specific pair of points under consideration will be 
denoted (l,j,k.“l) and (m,n,k). 

a. Check to see if any of the performance constraints 
are violated. These are: 

1) Maximum fuel flow for climb or cruise 

2) Descent prior to ascent, if 

and ALT (i,j,k-l) = 0 

3) Deceleration prior to acceleration if Vj > Vn 
and VEL(l,j,k-l) = 0 

4) Climb gradient ((Hi“Hni)/(D|^-Dj^.l) ) 

greater than 10%. This can easily be changed to 
a function of altitude. 

5) W(m,n,k) > 90,000 (indicates performance 
constraints violated to reach point (m,n,k>) 

b. If any of the above are violated: 

B(m,n) = 99,999 

c. Otherwise fuel burn from (i,j,k-l) to (m,n,k) is 
calculated and denoted B(m,n). 
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d. 


Thence 


W(i,j,k-1) = minm^n + W(m,n,k)]. 

If the optimum burn B(m,n) is less than 2% 
different from fuel burn with constant 
altitude/velocity, then latter is retained 
(i = m, j = n) . 

e. Flags are set; 

OH(i,j,k-l) = min(m) index 
OV(i,j,k-l) = min(n) index 
ALT(i,j,k-l) = 0 if > Hi 
VEL(i,j,k-l) = 0 if Vn > Vj 

f. Steps a to f are repeated for each (i,j,k-l). 

3. Decrement distance by one index and repeat step 2. If 
new Dj^-i = 0 (that is k -1 = 1 ), step 2 is only done 
for point ( 1 , 1 , 1 ). 

4. The optimum path is determined by path flags, iteratively 

a. Distance = D 2 

Density Altitude = H(0H(1,1,1)) 

Velocity = V(0V(1,1,1)) 

b. Distance = D 3 

Density Altitude = H[0H(0H(1,1,1) ,0V(1,1,1) ,2) ] 
Velocity = V[0V(0H(1,1,1)0H(1,1,1),2)] 

c. etc. 
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5. The density altitude is converted to pressure altitude 
before output by 

Hp = (TempH - 59. )/(“.003566) 

where : 

Tempjj = temperature at density altitude H 

Tempjj = C5 + C^h 

Hp =» pressure altitude in feet. 

The calculation steps are summarized in Figure 4-8. 
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Initialize 

Matrix 

Elements 


Calc. Burn From 

to (1,1, P) 
and calculate 
Weight (i,j,P-l) 


Set distance 
Index 


Calc. Burn (B(m,n)) from 
(i,j,k-l) to (m,n,k) if 
performance constraints aren’t 
violated. If k “ 2 then calculate 
burn only to point (1,1,1) 


Weight (i,j,k-l) = min^^^[B(m,ti) +W(m,n,k)] 
Set flags appropriate to indices m,n 


Decrement Distance Index 
k = k-1 


FIGURE 4-8 

MATRIX CALCULATION STEPS 










5. INFLIGHT UPDATE 


After the optimum flight path has been determined by the dynamic 
programming technique an option exists for a user Inflight 
update. If this option is selected, the fuel burn along the 
flight path, either indicated by the user or as determined by 
the dynamic program, is calculated. Also available is an option 
to update departure weight or aircraft weight at a distance node 
point . 

The steps through the inflight update are logically defined 
below (and in Figure 5-1) where OPH(k) and OPV(k) are the 
optimum altitude and velocity respectively at distance D(k) as 
determined by the dynamic program. 

1) Allow the user the option to set departure weight. 

Otherwise the weight at the departure point (W(l,l,l)) 
computed by the dynamic program will be used. 

2) Set the distance index k equal to 1 (departure point). 

3) Prompt user for distance node (L) where update is 
desired. (Note - only the weight can be changed at the 
departure point not the altitude or velocity). 

4) If the distance node for update (L) equals the node at 
the end of the flight segment (k+1) then the user will 
be prompted for the velocity and density altitude which 
will be set equal to OPV(k+l) and OPH(k+l) 
respectively. The user will also be offered the option 
to update the weight at the distance node (k+1). 

5) Compute the fuel burn from distance (k) to distance 
(k+1). 

6) If the user has updated the weight at this distance node, set 
the weight equal to the user input; otherwise, decrement the 
weight by the fuel burn of the segment, 

7) Increment distance index k = k+1. 

8) Steps 4 to 7 are repeated until either k equals the 
arrival point node or k=L. 

9) If k does not equal the arrival point node then repeat 
steps 3 to 8. 
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FIGURE 5-1 

UPDATE CALCULATIONS 
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6. CONCEPTS FOR EXPANSION 


The restricted search version of the dynamic programming 
approach, as described in Section 4, can be expanded to provide 
greater accuracy - 

Accuracy of the solution and speed of the computation can be 
enhanced by further consideration of factors such as those 
listed in Table 6-1. Depending on the description of the 
individual problem, the number of possible values for distance, 
altitude, and velocity can be different. An initial capability 
has been implemented in the current algorithm. The current 
algorithm can also be extended to consider unequal Increments in 
the division of distance, altitude, and velocity. This permits 
more efficient usage of fixed computational resources in those 
areas of distance, altitude, and velocity that are closest to 
the likely range of the optimal flight path. 

Perhaps the most critical improvement is in considering the 
limits of altitude and velocity as a function of distance. 
Clearly there are limits on climb and descent rates that can be 
reflected in the choice of altitude and velocity. Minimum 
values may be derived from pilot or file data on minimum en 
route altitudes and the like. Similarly, the number of 
candidate altitudes and velocities can be made a function of 
distance. 

Several additional features or algorithm capabilities could be 
considered. These could include 

o Time constraints 

o In-flight change of segments 

o Alternate airport planning 

o Point of no return indication 

o Adaptive aircraft constants to reflect 

- errors in data input (especially weather) 

- aircraft specific characteristics 

o Automatic Interfaces between ground computer and 

- weather sources 

- geographic data (to define location, minimum 
altitudes at standard navigational fixes) 


29 


I 



(jO 

O 


TABLE 6-1 

DESIGN VARIATIONS 


• For Fixed Computational Resources, Accuracy and Speed 
Solution Can Be Enhanced by Careful Consideration of 

— Vector Limits Different on Distance, Aititude, Veiocity 

— Unequai Distance, Aititude, Veiocity Segments 

— Calcuiation of Min/Max Altitude and Velocity as Function of Distance 



The results of the current algorithm should be checked and 
verified. For this purpose, it may be useful to develop a form 
of the algorithm for a mainframe computer, to provide for 
enhanced computational power and accuracy. In any case, the 
algorithm needs to be flight checked against actual aircraft in 
typical flight regimes* 

The current algorithm, based upon an Apple II computer, is not 
transportable to the cockpit. As such, the initial update 
capability described in Section 5 is not readily usable in a 
practical operational environment. To do this, another approach 
to a solution will be required. Section 7 describes an approach 
to a dynamic programming solution that restricts the moves that 
can be made in altitude and velocity between distance 
increments. This solution should be implementable on a 
hand-held programmable Calculator. 
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7. CONCEPT FOR IMPLEMENTATION 


Sections 4 and 5 defined the application of a restricted search 
dynamic programming technique to the problem of determining the 
unconstrained minimum fuel aircraft trajectory along a 
predetermined track or route. This section defines a 
modification to the technique that reduces the computational 
complexity and the computer memory requirements. The resulting 
technique then has the potential for implementation on a 
programmable calculator. 

7.1 Background 


The concept of utilizing a programmable calculator to perform 
pre-flight planning, inflight update and post-flight performance 
evaluation has been utilized in one case as part of the E3A 
AWACS fuel conservation and mission enhancement program. Figure 
7-1 depicts a HP-41C programmable calculator similar to that 
used in the application in conjunction with an appropriate 
overlay. The software represented mainly an attempt to minimize 
the amount of fuel burned on each trip segment in the presence 
of actual winds and temperatures and with a total trip time 
constant. This program Incrementally calculated each maximum 
range increment time and then adjusted this increment time so 
that the total time was consistent with a time constraint. It 
was recognized that, at best, this concept represented a first 
generation attempt at constrained trajectory optimization. The 
concept was flight tested and judged to be a success because it 
caused the aircraft to arrive at a geographic point at a precise 
time and with a fuel saving in the range of 5 to 7 percent. The 
next logical step was to upgrade the software so that it truly 
calculated an unconstrained minimum fuel trajectory. 

7.2 Algorithm Concept 


The algorithm described in Sections 4 and 5 must undergo a 
transformation in order to be compatible with the capability of 
the HP-41CV programmable calculator. A major simplification 
could be effected by the determination of relationships between 
altitude/distance and velocity/distance such that all feasible 
values could be examined in terms of a planar matrix. 

As given in Table 7-1, distance is used to determine the number, 
range, and value for the altitude and velocity points in the 
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'"'computes Between Way-' 
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Speed Command for 
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Mission Time 
Rec^uirements 
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Direction to 
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Static Air 
Temperature • *' 


•Clears Entry Error y’ 



Input Weather and 
Mission Track Data 
Output Min Fuel 
Mission Plan on 
Hard Copv and 
Mag Cards 


Computes Time to Leave 
Orbit With Fuel to Return 


CO 

CO 


FIGURE M 

PROGRAMMABLE CALCULATOR SAMPLE FORMAT 




TABLE 7-1 

TOWARD A REVISED ALGORITHM 
(RESTRICTION ON MATRIX MOVES) 


Using Altitude and Velocity as Function of Distance, Can Develop 
Revised Algorithm 

Distance Determines 

— Range of Altitude, Velocity Points 

— Number of Different Values To Be Used 

— Selection of Specific Altitude, Velocity Points 

Dynamic Program Steps Restricted 

— Arrive at Each Permitted Altitude, Velocity Point by One of Two Paths 
Only 



dynamic programming matrix. The dynamic program is thus 
modified from the "restricted search" approach by the addition 
of restrictions on the steps or moves in the dynamic programming 
solution. This "restricted move" approach is then the basis for 
a potential technique, implementable on a programmable 
calculator. 


This technique is illustrated with the aid of Figure 7~2. Part 
(c) shows a 10 X 10 matrix with the horizontal axis representing 
the variations in a kinetic energy index the vertical 

index representing the variation in a potential energy index 
(I ), and the diagonals representing the variation in a 
distance index (I^). The lower left vertex represents the 
departure point and the upper right the arrival point in a total 
path profile. Figure 7“2(a) depicts an altitude (h) surface 
that is a function of a potential energy index (Ip) and a 
distance index (Id)« This relationship can be expressed as, 

hp.d = Id) 

The ridge of the surface represents the altitude that can be 
achieved by the aircraft operating at maximum performance limits 
from departure to arrival. Thus an aircraft, at any distance 
index, can never reach an altitude greater than that which is 
located on the ridge of the surface. An empirical relationship 
that can be utilized to generate the surface is: 


P,d 


= C^e 




’t- 








In the above equation are constants and I^j and Ip as 
defined. Additionally, 

= f(trip maximum altitudes ceiling) 

(C 2 - C 4 , Cy) = f(aircraft performance constraints 

as derived from manufacturer's data) 

C 5 and Cg = f(total trip distance) 
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(a) 


Potential 
Energy 
Index (ip) 


Departure 


Kinetic Energy Index (i>/) 



1 


Distance Index (i^) 


19 


(c) 


FIGURE 7-2 

OPTIMAL FLIGHT PLANNING RESTRICTED MOVE TECHNIQUE 


In the example of Figure 7-2, ranges from 1 to 19. Given 
an then Ip ranges from: 

max (1, I^j - 9) to min (10, I^j) 

As an example, this equation means that at a distance from the 
departure point represented by I^j = 6, there are 6 values of 
Ip that can be selected, or there are 6 nodal points on the 
I^j = 6 diagonal matrix line. The selection of an Ip value 
defines a nodal point in the matrix and the equation then 
defines a unique altitude. This altitude relationship 
determines a family of ridges on which exists a possible optimal 
solution. An additional observation is that the maximum ridge 
level that can be reached in a segment is determined by the 
preceding part of the optimal solution. An explicit form for 
the lip^d function is given in Appendix E. 

Figure 7-2(b) shows a similar relationship to Figure 7-2(a) 
except that it concerns velocity (V), distance index (Ij) and 
a kinetic energy index (I^)* relationship is, 

V = fdd, 1 ^'). 

However, the maximum ridge generated does not represent a 
maximum velocity ridge; but, rather represents a maximum 
efficiency velocity ridge. The basic assumption is that the 
most fuel efficient velocity at an altitude is the best velocity 
to be selected if that velocity and altitude is a part of a 
total optimal solution. 

The equation has a similar form as that for the altitude case 
and is as follows: 

'"d,k = ^io® Li-^ -• 

In the above equation are constants, and I^j and Ij^ as 
defined. Additionally, 

Cio = f(trip maximum velocity) 

(Cii - C]^ 3 » = f (aircraft performance constants 

as derived from manufacturer’s data) 

and - f(total trip distance) 
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Thus, by means of the velocity and altitude functions It Is 
possible to compute an altitude and velocity at each nodal point 
after selecting an Ip, and Ij. This matrix Is pictured 
In Figure 7-2(c). 

7.3 Dynamic Program Calculations 

The first step Is to systematically compute a velocity and 
altitude at each nodal point In the matrix. The velocity at 
each nodal point is modified to ground velocity by the use of 
upper winds and temperatures data that have been entered. The 
altitude is with respect to air density and will be modified to 
pressure altitude only at time of output. 

The next step is to compute the total minimum fuel consumption 
required to go from any nodal point in the matrix to the 
departure point. Although the matrix presentation in Figure 
7“2(c) seems different than that of Section 4, many similarities 
exist . 

First fuel burns are calculated from the two points of = 18 
to the single arrival point (I^j = 19). This is a function as 
well of Hpjd> ^k,d» weight and aircraft constants. We then 
calculate the fuel burn (based on updated weight) from each of 
the three points at - 17 to the one or two points at 

= 18 to which they are connected. For each point in the 
matrix there are only one or two points at the previous 
(smaller) leading to the reference point. Similarly there 
are only one or two points leading from the reference point to 
the next (larger) I^j value. The observation is critical to 
the concept of a retricted move approach to the dynamic program. 

The process outlined above is repeated "backwards" through the 
matrix. At each step the flags to indicate the partial optimal 
path are stored as in the concept in Section 4. When the 
process concludes at = 1, the optimal path is retrievable 
from the flags. 

Some simplification in the calculations can be achieved by 
recognition of the various performance constraints as outlined 
in Section 4. 

A considerable savings in calculation is possible. For a sample 
of fdmax = 19 and = ^kmax = shown in Figure 

7-2, a maximum of 180 fuel burn calculations are required. For 
the approach of Section 4, a maximum of 160,200 fuel burn 
calculations could be required. (This number is drastically 
reduced with the consideration of performance constraints.) 
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It is emphasized that the technique has not been programmed and 
flight tested and thus should be considered preliminary until 
this work has been completed. 

7.4 Functional Description 


Section 4 describes the optimizing technique that computes the 
minimum fuel consumption path profile. Additional programs 
provide input and output data handling in terms of data storage 
organization. Since these programs are described in Sections 3 
and 5, they will not be repeated. Some details of the 
input /output process for the programmable calculator are 
contained in Appendix E. 
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APPENDIX A 


FUEL BURN MODEL 

A.l Turbojet Equation 

The basic equation for turbojet aircraft was derived in 
Reference 1. Much of this remains the same, and will not be 
repeated here in detail. 

A. 1.1 General 

We will use the following general definitions: 

W^ = fuel flow rate 

F = thrust 
n 

E^ = thrust work 

D = drag 

Ejj = drag work 

KE = kinetic energy 

PE = potential energy 

T = time in segment 

d = ground distance in segment 

= true velocity (average in segment) 

B = fuel burn 
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The general energy balance equation is 


Ex = Ep +APE +AKE 

We note also the following relationships for the case of no 

E_ D dx DV dt 

Do o T 

E, =/^ F dx * dF = TV„F 
Ton n T n 

B 3s Wf T 


A. 1.2 Original Derivation 

From Reference 1, the basic expression for Ep has been 
derived, and is 


E 


D 


pS TV 
1 1 ^ V 

2 


3 

T 


+ 


2R2K2W^T 


where 

R^ are flap/gear configuration drag multipliers 

are aircraft constants 
p is atmospheric density 

S is wing area 

W is aircraft weight. 


(A-l) 

wind: 

(A-2) 

(A-3) 

(A-4) 

(A-5) 
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A satisfactory form Is then found to estimate Fn/^f* This 
function is used since data is typically arrayed as: 




The empirical relation that defines this ratio is 


9 


(A-6) 


where 


are constants 
h is altitude - 


This yields then 


L / A_ 


Finally, we recognize 


= h (""xa - ^Ti) 


where g is gravitational acceleration. 
PE = WCh^ - h^). 


(A-7) 


(A-8) 


(A-9) 


We thence have 


B = W^T 


Ejj +AKE +APE 

e“/w^t 


(A-10) 


This, with the values of the preceding equations is the turbojet 
fuel burn equation. 
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A. 1.3 Low Altitude Enhancements 


The original solution was developed for climb/crulse/descent , 
for above 2000 feet and 250 Kn. Use of the fuel burn equation 
outside of these limits required a modification to reflect the 
higher thrust values for the low altitude /low speed regime. 
Thus, for thrust values in excess of Ky, we add to the fuel 
burn equation 


LAMl = KgFjj + Kq 


(A-11) 


where is derived from original fuel burn equation, noting 
it is of the form 

Erj, TV^Fn = Ej) +AKE +APE. (A-12) 

The detailed form of the equation is given in Figures A-1 and 
A-2. 


A. 1.4 Actual Winds and Temperature 


In the event actual winds and temperature are input into the 
fuel burn equation, account must be taken of the differences 
between ground and true velocity. 

True velocity (V;f) is used in describing the motion of the 
aircraft in the air mass, and determines the instantaneous fuel 
flow. Ground velocity (Vq) reflects the motion of the 
aircraft with respect to the ground, and determines aircraft 
progress along the route of flight. Ground velocity is defined 
as the magnitude of the ground velocity vector (Vq) defined in 
equation A-1 3. 

Vq = - WV (A-13) 

where Vq = vector with magnitude ground velocity and 
direction of true heading 

V-ji = vector with magnitude true velocity and 
direction of true course 

WV = vector with magnitude wind velocity and direction 
that the wind is blowing from. 
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B = 


+ T(LAMl) 


V^A 


where 


2 2 8tI’t2 

^ w T 


Vti1 + 


TV (*’2 ■ \) 


A = K ^ ^ + K 

A + K^2 


hj + 


+ K 


'*>1 + "2 


13 


+ K 


14 


LAMl 


if F < 
n 7 


KqF + otherwise 

on y 


R^jR^ contained in Figure A-2 


FIGURE A-1 

FUEL BURN EQUATION FOR TURBOJET (ISA CONDITIONS) 

Ui 



fGU^F^ + GU 2 F^ + GU^F + 1 (gear up) 

^ [gD^F^ + GD^F^ + GD^F + GD^ (gear down) 

= FDM^F^ + FDM^F^ + FDM^F + 1 
where 

F Is flap angle In degrees 

GU^, GDj^, and FDMj^ are aircraft constants 

FIGURE A-2 

FLAP AND GEAR DRAG MULTIPLIERS 
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In addition, actual conditions should be reflected in 
calculating density as a function of average temperature 
(T, in degrees Fahrenheit) and flight level (hp). 

Density below flight level 36089 is given by 

Pl, - 4 5 -17 67% T " • 0000068 (A-14) 
and above flight level 36089 it is given by 

Pj, = 45976^+ T ((36089 - hp)/53.35)] (A-15) 

where 

T = degrees Fahrenheit. 

The resultant fuel burn equation is given in Figure A-3. 

A. 2 Turboprop Equation 
A. 2.1 General 


Section A.l presented the equation form that can be utilized to 
compute the estimated fuel burn of a turbojet aircraft. Before 
developing a similar equation that will apply to propeller 
aircraft, it is relevant to discuss those differences in 
powerplants that account for equational differences. 

In the case of turbojet powered aircraft, the thrust force is 
generated reactively by the mass flow velocity of the jet nozzle 
gases in combination with the intake compressor by-pass 
airflow. The thrust is a function of the difference between the 
aircraft’s velocity through the air and the mass flow velocity 
through the engine nozzle. In turboprop powered aircraft the 
thrust force is derived mainly through the propeller which the 
turbine operates in addition to the compressor. The power 
generated is essentially constant with speed and the thrust 
decreases with speed. Thus the equational form of the fuel burn 
equations for the propeller aircraft and turbojet aircraft are 
different . 
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+ T(LAMl) 


B = 


TV^F^ 

VtA 


where 


F 

n 




2R2K2W 

Phv? 


w 

+ — (V 
gT T2 




A;LAM1; previously defined in Figures A-1 and A-2 
a is altitude dependent density 


T is derived as Ground Distance 
Ground Velocity 


FIGURE A-3 

FUEL BURN EQUATION FOR TURBOJET 
(ACTUAL WINDS AND TEMPERATURES) 



A. 2. 2 General Derivation 


It is convenient to concentrate on the relationships of power, 
since power relates to the shaft output of the engine. Data is 
readily available to establish a curve fit between P and Wf of 
the form: 


Wf = aP + be*^^ 


(A-16) 


where 

a,b,c are constants 
P is output shaft horsepower. 

Thrust and shaft horsepower are related by 


V^F 
T n 

np550 


(A-17) 


where r|p is the propeller efficiency, and is approximated by 
(.82), 

Substituting and letting 

Ki 5 = a/(Tip 550) 

Ki 6 = b 

%7 = 
results in 




Ki7h 


(A-18) 


We take the derived value of from equation (A-12) 


F = 
n 


+AKE +APE 
_ 


(A-19) 
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Finally 


h 

B = W^T = ^6*^^ * (A-20) 

The resultant fuel burn equation is given in Figure A-4. 

A. 2. 3 Actual Winds and Temperatures 


Actual winds and temperatures are reflected in the fuel burn 
equation analogously to the turbojet case. The resultant is 
given in Figure A-5. 

A. 3 Turbocharged Piston Equation 


A. 3.1 General Derivation 


As in the turboprop aircraft, the power is converted into thrust 
by the propeller. In a turbocharged engine, constant power 
output can be maintained to a certain altitude. 

Turbocharged piston engine aircraft normally use a lean fuel 
mixture during cruise and descent and richer fuel mixtures 
during takeoff and climb conditions. This can result in 
different fuel flows at the same altitude and velocity. When 
the lean fuel mixture is used, an empirical curve fit of fuel 
flow (W£l) Is a linear function of brake horsepower P: 


WfL = aP + b 


(A-21) 


From operating manuals it is possible to derive an empirical 
curve fit for fuel flow at richer fuel mixtures (Wf r) , 
assuming standard operating procedures. This function is 
quadratic instead of linear in P, to accommodate different fuel 
mixtures at different power settings: 


WfR = cp2 + sP + r 


(A-22) 


making the assumption that the richer fuel mixtures are only 
used when the aircraft is in climb. Equations (A-21) and (A-22) 
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B = K^5TV^F^ + K^^Te 


Kiyh 


where 


R,K^pS„V^ 2R,K^W^ ^ 

^n 2 _g ^2 gT ^^2 ■*■ TV ^**2 

^ w T 


previously defined 


FIGURE A-4 

FUEL BURN EQUATION FOR TURBOPROP (ISA CONDITIONS) 



Ln 

N3 


B = K^5TV^F^ + K^^Te 


Ki7h 


where 


F 

n 


RlK 


2R2K2W" 


+ — (V 
gT ^ T2 


T1 


)+|^^(h2-hi) 


previously defined 

p, is altitude dependent density 

T is derived as Ground Distance 
Ground Velocity 


FIGURE A-5 

FUEL BURN EQUATION FOR TURBOPROP 
(ACTUAL WINDS AND TEMPERATURE) 



can be combined using an exponential function, in the rate of 
climb (R/C), as a step function: 


=( 




](Wfi^) +(l - e 


K( 


<«fR> 


(A-23) 


With the constant K properly chosen, this exponential will be 
equal to 1 for climb, and 0 at all other times. 


Substituting 

R/C = V^(h2 - h^)/d 
P = V^F^/(np550) 
into equation (A-23) we have: 


W^ =6 


where 


(A-24) 

(A-25) 


/ aV_F ) 


2 2 

1 cv;;f 

sV„F \ 

1 T n 

I 4- b + (1 - 6 ) 1 

' T n 

+ ^ + rl 


^(np 550)2 

np550 j 


(A-26) 


, ^ [h - •’ih (**2 - 

Letting 

Ki5 - K 
*^16 “ 

Hi = ” 

Kig = c/(rip 550)2 
Ki 5 = s/(r,p 550) 
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and noting 

B * WfT (A-27) 

we have fuel burn as: 


» = ^ ^ ‘^17 > 
where defined In the 
The results are summarized 


+ T(l- 6 KK^g'^T^n + 
turbojet derivation 


+ '^ 20 > 

(A-28) 


in Figure A-6. 


A.3«2 Actual Winds and Temperature 

Accounting for actual winds and temperatures is analogous to the 
previous derivations. The resultant is presented in Figure A-7. 


54 



B = 6T(K^gV^F^ + + T(1 - 6)(K^gV2 


where 

^K^5(l/T)[|h2-h^[+(h2-hp] 


F 

n 


R^K 2R K W" 

-2^ P S/T+^-f-^iT 

pS vz 

^ w T 




FIGURE A-6 

FUEL BURN EQUATION FOR TURBOCHARGED PISTON 
(ISA CONDITIONS) 



Ln 

<y\ 


2 Jl 


B = 6T(K^gV^B^ + + T(1 F„ + 


20 ‘ 


where 


R K 2 

F = Pv,KK + 

n 2 ^h w T 


2R2K2W 

PhVt 


+ (V 
gT ^ T2 


'n> * Iv, <‘2- ‘i> 


6 previously defined 

previously defined 

is altitude dependent density 

T is derived as Ground Distance 
Ground Velocity 


FIGURE A-7 

FUEL BURN EQUATION FOR TURBOCHARGED PISTON 
(ACTUAL WINDS AND TEMPERATURE) 



APPENDIX B 


AIRCRAFT CONSTANTS 


B.l General 


This appendix contains aircraft constants for the fuel burn 
equation for both a turboprop aircraft (King Air 200) and a 
turbocharged piston aircraft (Cessna 421C). The use of these 
aircraft constants with the fuel burn equation and the techniques 
used to derive them are described in detail in References 1 and 
2, and Appendix A. Accuracy checks using the aircraft constants 
and cruise conditions are also given for both of these aircraft. 

B.2 Notation 


H = altitude (feet) 
V = velocity (true 
^%EF ” fuel burn, 
FBcal " fuel burn, 

^®PEF “ 


""CAL 


% = ^ — X 100 


FB 


REF 


airspeed in knots) 

reference value (pounds/hour) 

calculated value (pounds/hour) 


B.3 Beechcraft Super King Air 200 


Tables B-l and B-2 contain the aircraft constants and an 
accuracy check respectively for the Beechcraft Super King Air 
200 powered by two Pratt and Whitney PT6A-41 engines. The 
reference fuel burn values used in the accuracy check (Table 
B-2) were taken from the Super King Air 200 Pilot's Operating 
Handbook (Reference 5). 

B.4 Cessna 421C 


Tables B-3 and B-4 contain the aircraft constants and an 
accuracy check, respectively, for the Cessna 421C (Golden Eagle) 
powered by two Teledyne Continental Motors GTSI0-520-L engines. 
The reference fuel burn values used in the accuracy check (Table 
B-4) were taken from the Cessna 421 Golden Eagle Information 
Manual (Reference 6). 


57 


TABLE B-1 

BEECHCRAFT SUPER KING AIR 200 CONSTANTS 


= 303 


FDMi 

- 0 

= .025601A 


FDM 2 

= 0 

K 2 = .04241259 


FDM 3 

- -.0057 

K ^5 = .0000002692 


WFIDLE 

= .067 

Kjl 6 = .080443 



= 0 

= -.000034 


L 2 

- 0 

GDi ^ 0 


^3 

= 0 

o 

II 

CM 

Q 

O 


L 4 

= 0 

8 

w 

II 

o 


VNE 

= 289 

GD 4 = 2.3573 


VS 

= 75 

GU^ = 0 


MTOW 

= 12500 

GU 2 = 0 


OEW 

= 7755 

GU 3 = .01547 




Takeoff 

Climb 


Cruise 

0 

II 

< 

0 


0 

0 

II 

CM 

0 


0 

A 3 = .lOE-09 

-.44E-10 

- 

-.44E-10 

A 4 = -.42238E-05 

-.39419E-05 

- 

-.39419E-05 

A 5 = .28228 

.29681 

• 

29681 


Information from Beechcraft Super King Air 200 Pilot Operating 
Handbook and Airplane Flight Manual (Ref. 5) and Janets All the 
World^s Aircraft (Ref. 7) was used in deriving this set of constants 



TABLE B“2 

BEECHCRAFT SUPER KING AIR 200 ACCURACY CHECK 


Time = 3600 sec- 
Weight = 11,000 lbs- 


RPM 

H 

V 

Z?R.EF 

ISCAL 

A 

% 

1700 

0 

240 

946 

919 

-27 

-2-8 

1700 

2000 

24 5 

918 

903 

-15 

-1.6 

1700 

4000 

250 

892 

887 

- 5 

-0.6 

1700 

6000 

254 

866 

864 

- 2 

-0.2 

1700 

8000 

259 

846 

849 

3 

0.3 

1700 

10000 

264 

828 

833 

5 

0.6 

1700 

12000 

270 

812 

823 

11 

1.3 

1700 

14000 

275 

800 

807 

7 

0.9 

1700 

16000 

280 

792 

791 

-1 

-0.1 

1700 

18000 

280 

746 

747 

1 

0.1 

1700 

20000 

277 

690 

690 

0 

0.0 

1700 

22000 

275 

644 

642 

-2 

-0.3 

1700 

24000 

273 

598 

599 

1 

0.2 

1700 

26000 

269 

554 

552 

-2 

-0.4 

1700 

28000 

265 

514 

511 

-3 

-0.6 

1700 

29000 

263 

494 

492 

-2 

-0.4 

1700 

31000 

257 

458 

453 

-5 

-1.1 

1700 

33000 

250 

422 

417 

-5 

-1.1 

1800 

0 

245 

982 

956 

-26 

-2.6 

1800 

2000 

250 

954 

939 

-15 

-1.6 

1800 

4000 

255 

926 

923 

-3 

-0.3 

1800 

6000 

260 

902 

906 

4 

0.4 

1800 

8000 

265 

882 

889 

7 

0,8 

1800 

10000 

27 0 

866 

873 

7 

0.8 

1800 

12000 

2 76 

850 

863 

13 

1.5 

1800 

14000 

281 

838 

845 

7 

0.8 

1800 

16000 

283 

802 

810 

8 

1.0 

1800 

18000 

281 

746 

752 

8 

1.1 

1800 

20000 

279 

692 

700 

8 

1.1 

1800 

22000 

277 

646 

652 

6 

0.9 

1800 

24000 

274 

598 

603 

5 

0.8 

1800 

26000 

271 

554 

560 

6 

1.1 

1800 

28000 

267 

514 

518 

4 

0.8 

1800 

29000 

264 

494 

495 

1 

0.2 

1800 

31000 

258 

458 

455 

-3 

-0.6 

1800 

33000 

250 

422 

417 

5 

1.2 
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TABLE B-2 
(Concluded) 


RPM 

H 

V 

ZIref 

I?C AL 

A 

% 

1900 

0 

250 

1016 

995 

-21 

-2.1 

1900 

2000 

255 

988 

977 

-11 

-1.1 

1900 

4000 

260 

962 

960 

-2 

-0.2 

1900 

6000 

265 

940 

943 

3 

0.3 

1900 

8000 

270 

920 

925 

5 

0.5 

1900 

10000 

276 

904 

915 

9 

1.0 

1900 

12000 

281 

890 

897 

7 

0.8 

1900 

14000 

287 

876 

885 

9 

1.0 

1900 

16000 

286 

822 

828 

6 

0.7 

1900 

18000 

285 

770 

775 

5 

0.6 

1900 

20000 

283 

716 

721 

5 

0.7 

1900 

22000 

281 

668 

671 

3 

0.4 

1900 

24000 

279 

620 

625 

5 

0.8 

1900 

26000 

275 

574 

575 

1 

0.2 

1900 

28000 

271 

532 

531 

-1 

CM 

• 

O 

1 

1900 

29000 

269 

512 

511 

-1 

-0.2 

1900 

31000 

263 

474 

469 

-5 

-1.0 

1900 

33000 

256 

434 

431 

-3 

-0.7 


Mean = .06% 

Standard Deviation = .98% 
Variance = .95% 
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TABLE B-3 

CESSNA 421C (GOLDEN EAGLE) CONSTANTS 


S„ 

= 

215 


GUi 

= 0 



.0274935 


GU 2 

= .00041563 

K 2 

= 

.0415015 


GU 3 

= .030365 


= 

-.57564 


FDMj^ 

= 0 

He 

= 

2.636E-07 


FDM 2 

- 0 

H7 

= 

.0057675 


FDM 3 

= 0 

00 
1 — f 

= 

2.12978E-12 


WFIDLE 

= .01111 

Ki9 

= 

5.3845E-07 


1^1 

= 0 

^20 

= 

.0798467 


L 2 

= 0 

GD^ 

= 

0 


^^3 

= 0 

GD 2 

=S 

.0004145 


L 4 

= 0 

GD 3 

= 

.03042 


VNE 

= 258 

GD 4 

= 

1.9641 


VS 

= 74 

MTOW = 

7450 


OEW 

= 4426 



TAKEOFF 

CLIMB 


CRUISE 


= 

0 

0 


0 


= 

0 

0 


0 



-1.68E-10 

-1.68E-10 


-1.68E-10 

A 4 

= 

2.410128E-06 

2.410128E-06 


2.410128E-06 

A 5 

= 

.16033551 

.16033551 


.16033551 


Information from the Cessna 421 Golden Eagle Information Manual (Ref. 6) 
and Jane’s All the World’s Aircraft (Ref. 7) was used in deriving 
this set of constants. 
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TABLE B-4 

CESSNA 421C (GOLDEN EAGLE) ACCURACY CHECK 


RPM 

MP 

H 

V 

Z2ref 

ZiCAL 

A 

% 

1900 

32.5 

0 

186 

257 

254.3 

-2.7 

-1.0 

1900 

31.0 

0 

183 

246 

244.9 

-1.1 

-0.4 

1900 

29.0 

0 

177 

228 

227.1 

-0.9 

-0.4 

1900 

27.0 

0 

171 

211 

210.7 

-0.3 

-0.1 

1900 

25.0 

0 

164 

194 

193.1 

-0.9 

-0.5 

1900 

32.5 

5000 

195 

257 

255.7 

-1.3 

-0.5 

1900 

31.0 

5000 

191 

246 

244.0 

-2.0 

-0.8 

1900 

29.0 

5000 

185 

228 

227.4 

-0.6 

-0.02 

1900 

27.0 

5000 

178 

211 

209.6 

-1.4 

-0.7 

1900 

25.0 

5000 

171 

194 

193.3 

-0.7 

-0.03 

1900 

32.5 

10000 

204 

257 

255.4 

-1.6 

-0.6 

1900 

31.0 

10000 

200 

246 

244.5 

-1.5 

-0.6 

1900 

29.0 

10000 

194 

228 

229.0 

1.0 

0.4 

1900 

27.0 

10000 

187 

211 

212.3 

1.3 

0.6 

1900 

25.0 

10000 

179 

194 

194.9 

0.9 

0.5 

1900 

32.5 

15000 

214 

257 

256.1 

-0.9 

-0.3 

1900 

31.0 

15000 

210 

246 

245.9 

-0.1 

0.4 

1900 

29.0 

15000 

203 

228 

229.2 

1.2 

0.5 

1900 

27.0 

15000 

195 

211 

211.6 

0.6 

0.3 

1900 

25.0 

15000 

186 

194 

193.9 

-0.1 

-.05 

1900 

32.5 

20000 

224 

257 

255.2 

-1.8 

-0.7 

1900 

31.0 

20000 

220 

246 

245.9 

-0.1 

-.04 

1900 

29.0 

20000 

212 

228 

228.4 

0.4 

0.2 

1900 

27.0 

20000 

204 

211 

212.4 

1.4 

0.7 

1900 

25.0 

20000 

194 

194 

194.5 

0.5 

0.2 

1900 

32.5 

25000 

236 

257 

257.5 

0.5 

0.2 

1900 

31.0 

25000 

231 

246 

246.8 

0.8 

0.3 

1900 

29.0 

25000 

223 

228 

230.8 

2.8 

1.2 

1900 

27.0 

25000 

212 

211 

211 

0 

0 

1900 

25.0 

25000 

200 

194 

192.2 

-1.8 

-0.9 

1700 

32.5 

0 

177 

230 

227.1 

-2.9 

-1.2 

1700 

31.0 

0 

174 

219 

218.8 

-0.2 

-0.1 

1700 

29.0 

0 

168 

204 

202.9 

-1.1 

- .5 

1700 

27.0 

0 

161 

187 

186.0 

1.0 

.5 

1700 

25.0 

0 

154 

171 

170.7 

-0.3 

- .2 

1700 

23.0 

0 

146 

155 

155.2 

0.2 

.1 

1700 

32.5 

5000 

186 

230 

230.1 

0.1 

.04 

1700 

31.0 

5000 

181 

219 

217.0 

-2.0 

-0.9 
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TABLE B -4 
(Concluded) 


RPM 


H 

V 

Z2ref 

I®CAL 

A 

% 

1700 

29.0 

5000 

175 

204 

202.4 

- 1.6 

- 0.8 

1700 

27.0 

5000 

168 

187 

186.7 

- 0.3 

- 0.2 

1700 

25.0 

5000 

161 

171 

172.6 

1.6 

0.9 

1700 

23.0 

5000 

151 

155 

154.8 

- 0.2 

0.1 

1700 

32.5 

10000 

194 

230 

229 

- 1.0 

- 0.4 

1700 

31.0 

10000 

190 

219 

219.3 

0.3 

0.1 

1700 

29.0 

10000 

183 

204 

203.3 

- 0.7 

- 0.3 

1700 

27.0 

10000 

175 

187 

186.9 

- 0.1 

-.05 

1700 

25.0 

10000 

167 

171 

172.1 

1.1 

0.6 

1700 

23.0 

10000 

155 

155 

153.2 

- 1.8 

- 1.2 

1700 

32.5 

15000 

203 

230 

229.2 

- 0.8 

- 0.3 

1700 

31.0 

15000 

198 

219 

218 

- 1.0 

- 0,4 

1700 

29.0 

15000 

191 

204 

203.5 

- 0.5 

- 0.2 

1700 

27.0 

15000 

183 

187 

188.4 

1.4 

0,7 

1700 

25.0 

15000 

172 

171 

170.2 

- 0.8 

- 0.5 

1700 

32.5 

20000 

213 

230 

230.4 

0.4 

0.2 

1700 

31.0 

20000 

208 

219 

220.2 

1.2 

0.5 

1700 

29.0 

20000 

200 

204 

205 

1.0 

0.5 

1700 

27.0 

20000 

189 

187 

186.4 

- 0.6 

- 0.3 

1700 

25.0 

20000 

177 

171 

166.7 

- 4.3 

- 2.5 

1700 

27.0 

25000 

194 

187 

183.9 

- 3.1 

- 1.6 


Percent Differences: 

Mean = 0.16 
Standard Deviation = 0.64 
Variance = 0.40 
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APPENDIX C 


OPERATIONAL SCENARIO 


This appendix presents the concept of the operational scenario 
for the use of the flight planning aid. 

C.l Introduction 


The conservation of aircraft fuel starts with the Initial trip 
planning, and ends with a post flight evaluation of actual fuel 
usage. During preflight planning, the optimum fuel usage 
performance characteristics of an aircraft must be combined with 
route alternatives and upper air mass conditions to arrive at 
the flight path and operating parameters to achieve the maximum 
fuel economy. This Involves the evaluation of "trade off" 
situations such as operating at a less fuel efficient altitude 
in order to take advantage of an upper wind condition. 

This first level of fuel burn optimization is based upon 
pretakeoff route conditions. After becoming airborne, the 
actual weather conditions must be evaluated at reasonable time 
intervals and the required adjustments made in order to optimize 
the fuel conservation process. 

During the post flight review the actual versus predicted fuel 
usage can be evaluated. This is an important step because, 
during the flight, judgments were made based upon use of the 
technique. The review can be utilized to modify or adjust the 
individual technique so that future flights will benefit. 

Examination of these three steps indicates that the solutions, 
if performed manually, would *be cumbersome, time consuming, 
complex and error prone. Thus, the need for a flight planning 
aid points to some type of computer capability. The 
state-of-the-art of modestly-priced programmable calculators has 
reached the level that adequate computational capacity is 
readily available. 

The availability of relatively inexpensive minicomputers that 
are capable of interfacing with larger machines via modems 
provide a potential for performing complex track flight planning 
at virtually all airport locations. 
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C.2 Utilization Alternatives 


The operational utilization of the minimum peak consumption 
trajectory technique can be separated into four categories with 
each having different capabilities and data requirements. They 
are as follows: 

ALTERNATIVES 

PRE FLIGHT 

INFLIGHT 

I 

Programmable 

Calculator 

Programmable 

Calculator 

II 

Minicomputer/ 
Time Sharing 

Programmable 

Calculator 

III 

Minicomputer/ 
Time Sharing 

Performance 

Computer 

IV 

Performance 

Computer 

Performance 

Computer 


C.2.1 Alternative I 


The use of the programmable calculator for both preflight planning 
and inflight update and control represents the least cost 
implementation of the minimum fuel trajectory technique. During 
preflight planning (Table C~l) the pilot examines the winds and 
temperatures along various routes and chooses the route that, in his 
judgment, has the best conditions. The appropriate programs are 
loaded into the calculator by utilizing either bar code or magnetic 
cards. The total trip is divided into segments and the course, 
magnetic variation, and distance of each segment determined (Table 
C-2) . Next, the upper winds and temperatures at various altitudes 
for each segment are determined. The maximum acceptable operating 
altitude along with the landing weight (including alternate fuel 
reserves) is also determined. This data is entered interactively 
into the calculator. The output is a printed listing for each 
segment of distance, fuel consumption, speed, altitude, wind 
component, and temperature. Additionally, the complete matrix data 
is stored within the calculator and on magnetic cards for inflight 
use. 

On entering the cockpit for departure, the pilot loads his inflight 
programs and trip data cards. At departure the time and actual 
gross weight are entered. During climb the pilot can receive climb 
speed and power setting by entering altitude and temperature. At 
the first waypoint the pilot enters time, altitude, velocity. 
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a^ 

ON 


TABLE C-1 

PRELIMINARY FUNCTION CONCEPT 
STEP I: PILOT SUPPLIED DATA 


Route Of Flight Data 

• Determine Route Of Flight 

• Divide Into Approximately 20 Minute Flight Segments 

• Define Distance And Course For Each Segment 

En Route Weather 

• Define Wind Velocity, Heading And Temperature 

— Maximum Operating Altitude 
(Based Upon Aircraft Limitations, 

Icing, Turbulence, Oxygen Availability) 

Calculator Programs 

• Have Available Program Modules/Magnetic Tape With Particular 
Aircraft Constants Stored 

Or 

• General Program Module/Magnetic Tape/Paper Copy With A List 
Of Particular Aircraft Constants 



TABLE C-2 

PRELIMINARY FUNCTIONAL CONCEPT 
STEP II: PRE-FLIGHT PLANNING 


Initial Data 

• Aircraft ID, Engines 

• Empty Weight 

• Altitude Constraint 

Segment Data Input 

• Distance 

• Altitude 

• Temperature 

• Wind 

• Course 

Output Best And Next Best 

• Climb, Cruise, Descent 

— lAS/Mach 
— Power Settings 
— Altitudes 



temperature, and gross weight or fuel on board (Table C-3) . 
Also, at this point a projected minimum fuel arrival time will 
be displayed. The pilot can either accept or modify the 
value. The program then reprojects, incrementally, the 
remaining portion of the flight and displays a speed, altitude, 
and power setting command. 

The entry is reported at each way point throughout the flight. 
The time is continually updated and can be displayed upon 
command to the pilot. Additionally, a start descent time is 
calculated and displayed. At start of descent the pilot enters 
the time, altitude, temperature, and fuel on board at which the 
descent speed is displayed with the appropriate power setting. 
As in climb, the speed and power setting is displayed at entry 
of a time, altitude, and temperature. Upon landing, the time 
and fuel on board is entered. The data contained in the 
calculator is recorded on magnetic cards. 

The post flight analysis (Table C-4) is performed by the pilot 
entering the appropriate programs in combination with the 
flight data that was recorded on the data cards. A printout 
records the actual versus projected plans and allows the pilot 
to evaluate the trip fuel efficiency. 

C.2.2 Alternative II 


This alternative assumes that the ground or preflight planning 
is performed by a minicomputer that is linked to a VOR and a 
winds and temperatures data base (Table C-5). It also assumes 
that the computer has the capability to evaluate multiple 
tracks either iteratively or by utilizing a dynamic programming 
data base. The pilot enters landing gross weight, arrival 
point (VOR), maximum acceptable altitude, departure point, and 
aircraft /engine types. The system then provides a printed 
output as described in Alternative I, including a bar code for 
loading into a programmable calculator. 

The inflight and postflight operations are the same as 
described in Alternative I. 

C.2.3 Alternatives III and IV 


These alternatives are listed for the purpose of highlighting 
the point that the techniques described can be Implemented on 
any computer, providing that the required interface capability 
has been provided. 
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TABLE C-3 

PRELIMINARY FUNCTIONAL CONCEPT 
STEP III: IN-FLIGHT UPDATE AND CONTROL 


Updated Inputs For Way Points 

• Time 

• Altitude 

• Fuel Burn/Flow 

• Temperature 

Revised Output 

• lAS/Mach 

• Power Settings 

• Altitude 



TABLE C-4 

PRELIMINARY FUNCTIONAL CONCEPT 
STEP IV: POST-FLIGHT REVIEW INPUT AND OUTPUT 


To Be Determined Based On Pilot Needs 
Objectives 

• Compare Best To Actual Fuel Conservation Performance 

— stored In Calculator 
— Reflects Encountered Conditions 

• Present Results So That Future Performance Will Be Improved 



TABLE C-5 

OPERATIONAL SCENARIO 


Prefiight 

— Time Share or Separate Computer at Fixed Base Operator 

— (Later) Possible Automated Weather Input 

— Initial Optimal Solution 

Inflight 

— Programmable Calculator, Initialized With Ground Solution 

— Inflight Updated Inputs and Solution 



jlA U 


SAMPLE PROBLEM 


This appendix contains a sample optimization problem done on 
the Apple II computer. This problem consisted of determining 
the optimum path profile for a Super King Air 200 aircraft 
during a 360 nml flight. 

D.l Inputs 


The actual Inputs for this problem are given in Table D-1 (as 
echoed out to the printer as output). The altitude modes were 
calculated by the program while the velocity modes were 
designated by the user. 

D.2 Matrix Definition 


The matrix nodal points in altitude, velocity, and distance are 
defined by either the user or the program. In the sample 
problem the velocity nodes are defined by the user while the 
altitude nodes are defined by the program. Since the number of 
distance nodes is greater than the number of waypoints, the 
program embeds distance nodes between waypoints. The altitude, 
velocity, and distance nodes defined for the sample problem 
were printed out and are contained among the output figures in 
Table D-1. 

D.3 Optimal Path Profile 


Once an optimum path profile has been determined, the path is 
printed in the form shown in Table D-1, along with the 
departure and arrival weights. The altitudes printed out 
differ from the altitude nodal points since the nodal points 
are density altitudes while pressure altitude is given. The 
program also prints a list of flight segments using the optimum 
path containing start and end distances, velocities (TAS), and 
pressure altitudes for each flight segment, and the fuel burn 
over the segment as shown for the sample problem in Table D-1. 

D.4 Outputs 


The outputs for this sample problem are shown in Table D-1 and 
have been described previously. 
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TABLE D-1 
SAMPLE OUTPUT 


FLIGHT TITLE 
? SftMPLE FLIGHT 

INPUT SER - CEILING i-DEP -RLT . -.DEP . VEL ?RF:R -RLT ,RRR ,VEL 
? 33000- j5000, -1135. ^5000 - 135 . 

NO - OF RLT - NODES' -.VEL -NODES^DIST - NODES jWhVPDINTS? 

? 1 0)^1 0^1 0)i5 

INPUT RLT I TUBE NODES •Si' -'-N) 

? N 

INPUT VELOCITY NODES-:. Y.-N> 

? Y 

INPUT VELOCITY NODES 

? 1:35- .*155- .*175- .*190, j£00- .*210- .*22 0- j230- .*250- .*27 0- 
LRNDIHG GROSS WEIGHT 
7 11250- 


WRYPOINT NUM.BER 1 

DISTRhCE= 0.000000 
HERDTNG= 40-0000 


RRIRTIDN= 

:3 - 0 0 ij 0 




RLT=15 000,0 

WD= 90.0 

UV=3 0,000 

TEMP= 

5 , 5 0 

RLT=25000 , 0 

WD= 70-0 

l.i .1 V =■ 4 ij , U U U 

TEMF- 

-30. 2 0 

RLT=35000 . 0 

WD= 5C.0 

!,jV=7 0 . 00 0 

TEMF- 

-65 .80 


WRYPOINT NUNPER 

D I ST RN C E = 5 4 - 0 1'.i 0 0 0 0 
HERDIN6= 120,0000 
VRRIRTION= 3. 00 CO 


RLT=15 00 0 . 0 

WD-=16 0.0 

W V=55 . 0 0 0 

TENF = 

2 - 0 0 

RLT=25 0 00 , 0 

WD=13C,0 

WV=60 , 000 

TENP= 

-22 . 0 0 

RLT=35 00 0 , 0 

WD=c 3C - 0 

WV=7 0.000 

TEMF- 

-58 . 00 


WRYPOINT NUMBER 

“I 




DISTRNCE=115, 

0 0 0 0 0 c 




HERDING= 

8 0 . 0 0 0 0 




VRRIRTIDN= 

5 . 0 0 0 0 




RLT=15000 , 0 

WD- 1 3 0 , 0 

WV=40 , 000 

TEMP = 

-1,0 0 

RLT=25000 . 0 

WD=18C,0 

WV=50 , 000 

TEMP= 

-25 . 00 

RLT=35000.0 

WD=21G-0 

WV=6 0.000 

TEMF== 

-59 .00 
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TABLE D-I 
CContinued) 


WRYPDINT NUMBER 4 

D I STRNCE=c:65 - 0 0 0 0 0 u 
HERDIM6= 90.0000 


VRRIflTIDH= 

8.00C0 




RLT=15000.0 

WD=e?G„D 

WV=3!:i.Ci00 

TEMF= 

-3 . 00 

RLT=£5000.0 

UD=cb C . 0 

WV=4 0 .000 

TEMP= 

-36 . 0 0 

RLT=35000.C 

WD=c9 0 - 0 

WV=50.000 

TEMF= 

-63 . 00 


WRYPDIMT MUMPER 5 

DISTRHCE=315.0000C0 
HERDIMG= 90.0000 

VRRIRTIDM= S.OOro 

RLT= 1 5 0 0 0.0 Wri=3 0 0 . 0 

WV=10.000 

TEMP= 

5.50 

RLT=£5 0 0 0 . 0 WD=£9 C , 0 

WV=£ 0.000 

TEMF= 

-30.30 

RLT=55 0 0 0 . 0 WD=87 0 . 0 

WV=40.000 

TEMF = 

-65 - 3 0 


RLTITUHE 

MODE 

1 = 

5 Cl 0 0.0 0 0 0 

RLTITUDE 

MODE 

o— 

u, — 

Sill .1111 

RLTITUDE 

HDDE 

3= 

1 1233 - £333 

RLTITUDE 

MODE 

4= 

14333.3333 

RLTITUDE 

MODE 

5= 

17444 .4444 

RLTITUDE 

MODE 

6*= 

30555.5556 

RLTITUDE 

MODE 

■J7 — 

i£-^666 -6667 

RLTITUDE 

HDDE 

3= 

36777.7773 

RLTITUDE 

HDDE 

9= 

39SSS .S8S9 

RLTITUDE 

HDDE 

10= 

3:5 0 0 0 . 0 0 0 0 


VELOCITY 

HDDE 

1 = 

135.0000 

VELDCITY 

MODE 

3= 

155.0000 

VELOCITY 

HDDE 

3= 

175.0000 

VELDCITY 

MODE 

4= 

190.0000 

VELDCITY 

HDDE 

5= 

300.0000 

VELDCITY 

MODE 

6= 

310.0000 

VELDCITY 

HDDE 

7= 

330.0000 

VELDCITY 

HDDE 

8= 

33 0.00 0 0 

VELDCITY 

HDDE 

9= 

350.0000 

VELDCITY 

HDDE 

10= 

370.0000 
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TABLE Dr^l 
CContinued) 


DISTANCE NODE 
DISTANCE NODE 
DISTANCE NDDE 
DISTANCE NODE 
DISTANCE NODE 
DISTANCE NODE 
DISTANCE NODE 
DISTANCE NODE 
DISTANCE NODE 
DISTANCE NODE 


1= 

0.0000 

2= 

£7.0000 

3= 

54.0000 

4= 

33.5000 

5* 

113.0000 

G= 

151 .0000 

r' “ 

139.0000 

3= 

££7.0000 

9= 

£65.0000 

10= 

315.0000 


ARRIVAL WEI6HT= 11 £5 0.0000 
DEPARTURE WEIGHT* 117S4.3645 


□PTIMUM PATH PROFILE 


NODE 

DIS. 

ALT . 

v*;:ktrs 

1 

0 . 0 0 

5010. £3 

1 35 . 0 0 

£ 

£7 . 0 0 

£03£9.35 

1 9 0 . 0 0 


54.0 0 

£7949-15 

1 9 0 . 0 0 

4 

S3 . 5 0 

3 0731 .45 

£00 . 00 

5 

113. 0 0 

30996.45 

£00. 0 0 

6 

151 .00 

31164.70 

£ 0 0 . 0 0 

7 

139.00 

3 1 33£ . 96 

£ 0 0 . 0 0 

3 

££7.00 

31-501 .££ 

£00.00 

9 

£65 . 0 0 

£3943.40 

£70.00 

10 

315.00 

5010. £3 

135.00 
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TABLE D-1 
(Concluded) 


□FT I MUM PftTK FUEL BURN 


D1 

D£ 

VI 

V 

0.00 

£7.00 

135.00 

190 

£7.00 

54-00 

190.00 

190 

54.00 

S3. 50 

190.00 

£00 

S3. 50 

113.00 

200.00 

£00 

113.00 

151 .on 

£00.00 

£00 

151 -00 

139.00 

£00. 0 0 

£00 

139.00 

££7.0 0 

2 0 0 . 0 0 

£00 

££7.00 

£65 . 00 

2 0 0 . 0 0 

£70 

£65.00 

315.00 

£70.00 

135 


HI HE B 


00 

5010 .es 

£0329.35 

109-62 

00 

£0329.35 

£7949.15 

72.91 

00 

£7949.15 

30781 -45 

51 .76 

00 

307S1 .45 

30996 -45 

40.52 

00 

30996 .45 

31164.70 

52-56 

00 

31164.70 

31332-96 

53.77 

0 0 

3 1332.96 

31501 .££ 

55 . 04 

0 0 

31501 .££ 

£3948.40 

40 . 02 

00 

2394S.40 

5010.2S 

6 0 . 0 1 



APPENDIX E 


PROGRAMMING DETAILS 


This appendix describes some of the programming detail 
necessary to Implement the concepts presented In the paper. 

E.l Apple II Algorithm 


This section describes programming details to Implement on an 
Apple II computer (In FORTRAN) the dynamic programming 
technique described previously. 

E.1.1 Storage 


Since the Apple II computer has limited storage, it is 
necessary to combine several arrays together to permit maximum 
usage of storage. In the Apple II program the OH and OV arrays 
described In Section 4.3 and the ALT and VEL arrays also 
described In Section 4.3 have been combined respectively. The 
OH and OV arrays are combined into an array lOP such that 

lOP (l,j,k) = OH (l,j,k)*100 + OV (i,j,k). 

The ALT and VEL array have been combined such that 

IFLAG (l,j,k) = ALT (1, j,k)*10 +VEL (i,j,k) 

E.l. 2 Input Data 


The following inputs to the Apple II will be supplied by the 
user. 

1. Flight Plan Title (25 characters or less) 

2. Effective Service Ceiling (feet) (Program will 
default to Aircraft Service Ceiling if exceeded) 

3. Departure Altitude (feet) (density altitude) 

4. Departure Velocity (knots) 

* 

5. Arrival Altitude (feet) (density altitude) 

6. Arrival Velocity (knots) 
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7. Number of Altitude Nodes (If the user wishes to input 
the altitude node points, then he ii^uts a ”Y" when 
prompted and then inputs the density altitude node 
points. These node points must be such that 

8. Number of velocity nodes (User can input velocity 
nodes in similar fashion to inputted altitude nodes.) 

9. Number of Distance Nodes 

10. Number of Waypoints (Must be equal to or less than 
number of distance nodes.) 

11. Landing Gross Weight (Must be greater than the 
aircraft operating empty weight else program will 
reprompt for Landing Gross Weight.) 

The program will then indicate the waypoint number and the user 
will then supply the following for that waypoint: 

12. Course Heading (degrees from North) 

13. Magnetic Variation (degrees from North) 

14. Distance to Next Waypoint (nautical miles) 

15. Number of altitudes at this waypoint for which 
weather data will be supplied 

16. Altitude (density altitude in feet) 

17. Wind Direction at Altitude (degrees from North) 

18. Wind Velocity at Altitude (knots) 

19. Temperature at Altitude (in degrees Farenheit) 

Once the optimization has been completed and the optimal path 
profile has been printed out the user will be asked if he 
wishes to update. If the user responds with a "Y” then the 
following inputs will be required: 

1. Departure Weight if the user responds in the 
affirmative to the prompt concerning updating 
departure weight. 
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2. Node Point where weight altitude, or velocity are to 
be updated. (This node point must be greater than 
previous node point updated or else the program will 
assume no additional updates will occur.) 

3. Updated Weight at node point if the user responds in 
the affirmative 

4. Density Altitude (feet) at node point 

5. Velocity (knots) at node point 
E.1.3 Outputs 

The following outputs from the Apple II computer to the user 
will normally Include: 

1. Flight Plan Title 

2. Effective Service Ceiling (feet in density altitude) 

3. Departure Altitude (feet in density altitude) 

4. Departure Velocity (knots) 

5. Arrival Altitude (feet in density altitude) 

6. Arrival Velocity (knots) 

7. For each user waypoint 

a. Waypoint Number 

b. Distance from Departure Point (nautical miles) 

c. Course Heading (degrees from North) 

d. Magnetic Variation (degrees from North) 

e. For each Altitude 

a) Altitude (feet in density altitude) 

b) Wind Direction (degrees from North) 

c) Wind Velocity (knots) 

8. Arrival Weight (W(1,1,P)) (pounds) 
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9. Departure Weight (W(l,l,l)) (pounds) 

(If the aircraft’s Maximum Takeoff Weight is exceeded 
a warning will be printed.) 

10. Optimum Path Profile 

This will include for each node point 

a. Node point (number to be used for updating) 

b. Distance from departure point (nautical miles) 

c. Pressure altitude (feet) 

d. Velocity (true airspeed in knots) 

11. Optimum Path Fuel Burn 

Using the Departure Weight (W(l,l,l)) the fuel burn 
will be computed along the optimum path. Outputs 
will include for each segment the start and end 
distance (D;^,D 2 ), the start and end velocities 
(V1,V2), the start and end pressure altitude (H1,H2), 
and the fuel burn (B) over the segment. 

12. Updated Path Fuel Burn 

If Update has occurred the fuel burn along the 
updated profile will be computed using either the 
optimum departure weight (W(l,l,l)) or the user’s 
departure weight; Outputs will include the same as 
with the Optimum Path Fuel Burn. In addition, both 
the departure weight and updated arrival weight will 
be printed. 

A sample output was included in Appendix D. 

E.2 Programmable Calculator Input Proce s s ing 


The preflight planning program for the progra mma ble calculator 
version has been divided into three subprograms or modules. 
These are: 

o Data Entry and Initialization 

o Flight Path Optimization 

o Data Output and Write Storage 



Each of these modules accepts data, performs computations, and 
outputs data either Internally or externally. The Flight Path 
Optimization module is common to the pre-flight; in-flight 
update; and, the post-flight evaluation. Of the three modules 
only the Data Entry and Initialization programming has been 
completed. Therefore, this section will deal only with this 
module. 

The purpose of the Data Entry and Initialization subprogram is 
to load all programs, aircraft specific constants, general 
constants, optimization algorithm specific constants, and the 
initialization program. Figure E-1 depicts a flow chart that 
defines the program flow along with explanatory comments. 

E.3 Programmable Calculator Optimization Constants 


The purpose of this section is to define the data that is used 
from the pilot’s handbook and to show how the h = f(I(j,Ip) 
and V = f(I^j,I]j^) equations' constants are determined. 

E.3.1 Pilot’s Handbook Data 


The following table depicts the data that was derived from the 
Beechcraft Super King Air 200 Pilot’s Handbook: (A 300 mile 

trip is assumed as an example; however, the non-trip specific 
constants will remain the same.) 


Distance 

(nmi) 

Altitude 

(ft) 

Velocity 

(knots) 

Id 

2.3 

1,000 

125 

0.146 

7.3 

10,000 

143 

0.462 

22.4 

22,000 

166 

1.42 

106.5 

35,000 

230 

6.74 

151.2 

35,000 

230 

9.57 

222.79 

22,000 

169 

14.10 

276.28 

10,000 

138 

17.49 

300.00 

1,000 

119 

19.00 


The above data describes a maximum performance climb, cruise, 
and descent based on the handbook data. 

E.3. 2 Constant Determination 


The constants contained in the altitude (h) and velocity (v) 
relationships are both aircraft type and trip dependent. In 
this typical example, a trip distance of 300 nautical miles 
will be used. 
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LOAD 

& 

START 


Subprogram loaded from magnetic cards 
which loads data registers with aircraft 
specific and other constants in addition 
to the initialization program. 


ENTER TRIP 
SPECIFIC DATA 


STORE DATA 


P + 1 


i + 1 


ENTER 
WAYPOINT 
i Data 


ENTER WAY POINT 
i, ALTITUDE 
P DATA 


COMPUTE WIND 
COMPONENT AT 
P ALTITUDE, 
WAY POINT i 



■ 




COMPUTE WAY POINT 
i WCP FUNCTION 
CONSTANTS 





COMPUTE WAY POINT 
i TEMP. FUNCTION 
CONSTANTS 



TO OPTIMIZATION 
PROGRAM 


- AIRCRAFT S ENGINE TYPE. 

- TOTAL TRIP LENGTH. 

- NUMBER OF WAY POINTS. (W WAYPOINTS) 

- EFFECTIVE AIRCRAFT CEILING. 

- LANDING WEIGHT (INCLUDING RESERVES). 

- DEPARTURE POINT ALTITUDE. 

- ARRIVAL POINT ALTITUDE. 


- DISTANCE TO NEXT WAY POINT 

- MAGNETIC COURSE TO NEXT WAY POINT 

- MAGNETIC VARIATION (EAST IS (-)) 

- i = 1 IS DEPARTURE POINT 

- NUMBER OF DATA ALTITUDES (J ALTITUDES) 


- WIND VELOCITY (KTS) AT P ALTITUDE 

- WIND DIRECTION (TRUE) AT P ALTITUDE 

- TEMPERATURE (°F)AT P ALTITUDE 


- WCP = -j^WIND velocity! C0S(C0URSE-WIND MAGNETIC) 

WIND MAGNETIC = WIND DIRECTION + VARIATION 
NOTE: 

- WCP = HEADWIND 


WCP = (CONSTANT 1) + (CONSTANT 2) (ALTITUDE) 

Utilize linear regression to 
compute constants 1 & 2 


TEMP(°K) = (CONSTANT 3) + (CONSTANT 4)(ALTITUDE) 

UTILIZE LINEAR REGRESSION TO 
COMPUTE CONSTANTS 3 & 4 


FIGURE E-1 

DATA ENTRY AND INITIALIZATION 
















Determination 


The relationship of the distance index (I^) to the ground 
trip distance (dg) is: 

dg 

where, is a constant for a specific trip, 

= 0.0633dg, 

for 19 points, and a dg, total = 300 nmi. 

Altitude Function 


As described in Section 7, the altitude (h) function is, 

hpd = f(Ip)f(I^j) (assuming separability) 

An equation form that will generate a ridge line in altitude as 
a function of is. 


h^ = ^^^d^ ” ceiling altitude 



-Cnid 

In the above equation the 1 - e term accounts for the 

maximum performance climb to the effective service ceiling 

-Cn |ld - 19[^n 

altitude. The term, e , represents a maximum 

range descent at minimum power condition. The other part of 


C (I “C ) 
n^ p n^ 


This 


the h ^ function is, hp = f(Ip) = e 
equation provides the cross section shape or variation in 
altitude with Ip. 


For the Super King Air 200, the constants have 
and are as follows: 

-0.1755(1 -5.5)r -0.69751 

h = 35,000 P 1 - e - e 


been determined 


-0.0175 |I,-19| 
d 


This example serves to demonstrate this method that is utilized 
to determine the h and v functions. Since the v function is 
similar, an example will not be presented. 
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APPENDIX F 


Symbol 

a 

ALT 

b 

B 

c 

Ci 

d 

di 

D 

Ed 

Et 

F 

FDMi 

En 

8 

GD 

GU^ 

h 

hp 

H 

Hp 


NOMENCLATURE 

Meaning 

Temporary constant 
Altitude (program variable) 
Temporary constant 
Fuel burn 

Temporary constant 
Constants 
Ground distance 
Distance steps 
Drag 

Drag work 
Thrust work 
Degrees of flap 
Flap/gear constants 
Thrust 

Gravitational acceleration 
Flap/gear constants 
Flap/gear constants 
Altitude 
Flight level 

Altitude (program variable) 
Pressure altitude 
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hp,d 

i 

Id 

Ik 

I OP 
IFLAG 

j 

k 

K 

Ki 

KE 

L 

LAMl 

M 

N 

OH 

OPH 

OPV 

OV 

P 

PE 


Ri 

R/C 


Altitude (Indexed by potential energy and distance) 

Density altitude index 

Index of distance 

Index of kinetic energy 

Index of potential energy 

Optimum index (program variable) 

Optimum flag (program variable) 

Velocity, index 
Distance index 
Constant 

Aircraft specific constants 

Kinetic energy 

Iteration counter 

Low altitude term 

Limit of density altitude index 

Limit of velocity index 

Altitude index 

Update altitude index 

Update velocity index 

Velocity index 

Output shaft horsepower (or limit of distance index) 
Potential energy 
Temporary constant 

Flap/gear configuration drag multipliers 
Rate of climb 
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s 


Temporary constant 


T 

Temph 

V 

Vg 

% 

Vx 

vi 

VEL 

Vd,k 

w 

Wf 

WfL 

WfR 

wc 

WE 

WN 

WV 

6 

P 

Ph 


Wing area 

Time (or Temperature) 

Altitude dependent temperature 

Velocity 

Ground velocity 

Vector with direction of true heading and magnitude 
ground velocity 

True velocity 

Vector with direction of true course and magnitude 
true velocity 

Velocity (program variable) 

Velocity (indexed by distance and kinetic energy) 

Aircraft weight 

Fuel flow rate 

Fuel flow rate (lean) 

Fuel flow rate (rich) 

Head wind component 
East wind component 
North wind component 

Vector with direction that the wind is blowing from 
and magnitude of wind velocity 

Ratio constant (in fuel flow) 

Propeller efficiency 

Atmospheric density 

Altitude dependent density 
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